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We present a measurement of the topological susceptibility in SU (3) Yang-Mills theory through the deconfine- 
ment phase transition. An improved operator is used for the topological charge density. A drop by an order of 
magnitude is observed from the confined to the deconfined phase. 



1. INTRODUCTION 

The topological susceptibility is defined in con- 
tinuum QCD as 



X= / d^x{0\T{Q{x)Q{0))\0), 



where 



Q{x) 
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Q{x) is the topological charge density and is pro- 
portional to the anomaly of the Ua{^) current 



d^j^ix) = 2NfQ{x). 



(3) 



The value of x for quenched QCD can be related 
to the masses of the pseudoscalar mesons, by use 
of expansion 
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X has been determined on the lattice ^ and 
is consistent with the prediction of eq. (^). This 
solves the so-called C^a(1) problem ||]. 

It is interesting to determine the behaviour of 
X at T 7^ , in particular across the deconfining 
phase transition, which in quenched QCD takes 
place at Tc ^ 260 MeV |^ . A general expectation 
is that the topological susceptibility should drop 
at Tc , since Debye screening inhibits tunneling 
between states of different chirality and damps 
the density of instantons. 

Previous attempts to study x through Tc by 
the field theoretical method were plagued with 
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large statistical uncertainties above Tc Also 
the attempts to use the cooling technique encoun- 
tered difficulties at finite T §,§. The idea of 
the cooling method is to measure the topo- 
logical charge on configurations cooled by locally 
minimizing the action. At T = and below Tc 
the method works well and instantons show up as 
plateaux of the topological charge versus cooling. 
Around Tc instantons become unstable, the plat- 
eaux disappear and the method becomes ambigu- 
ous §. 

In the present paper we drastically reduce 
the statistical fluctuations of the field theoreti- 
cal method determination, by using an improved 
topological charge operator as proposed in ref. 

We redetermine x at T = and study its be- 
haviour through Tc- The new determination of x 
at T = is consistent with previous determina- 
tions |l|-[l3). At finite T we obtain a good deter- 
mination of X- Our main result is that x drops 
at Tc by more than one order of magnitude. 

2. THE METHOD 

The field theoretical method is a straightfor- 
ward application of the basic rules and concepts 
of quantum field theory. 

At first a lattice discretization Ql{x) of the 
continuum topological charge density operator 
Q{x) is defined. Ql{x) is not unique: infin- 
itely many operators can be defined on the lattice 
which have the same continuum limit. 

In going to the continuum limit a proper renor- 
malization must be performed, like in any other 
regularization scheme. In quenched QCD Ql 
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renormalizes multiplicatively [|^. In formulae 

QL = Z{P)Qa'' + 0{a^). (5) 

As usual, (3 = 6/(7q. The topological susceptibil- 
ity can be defined on the lattice as 

XL^{Y.QL{x)QLm. (6) 

X 

The standard rules of renormalization then give 

XL = Z{pfa\ + M{p)^0{a^), (7) 

where M{(3) is an additive renormalization con- 
taining mixings of xl to other operators with the 
same quantum numbers and lower or equal di- 
mensions |ll5| . 

Z and M can be computed non-perturbatively 
on the lattice by the heating method ||l^ , |l6| , p7| . 
Renormalizations are induced by quantum fluctu- 
ations at the scale of the UV cutoff (lattice spac- 
ing scale), so they can be estimated by measuring 
topological quantities on an ensemble of configu- 
rations where the topological content is known 
and short-range fluctuations are thermalized. In 
order to do this one can start with a classical con- 
figuration on the lattice and perform few local 
updating steps on it until the short-range fluc- 
tuations are thermalized: by repeating this pro- 
cedure several times an ensemble of trajectories 
can be constructed, all starting from the same 
initial configuration. Starting from a large in- 
stanton and measuring on the ensemble, a 
plateau at the value of Z{(3) is reached with re- 
spect to the number of heating step performed. 
Similarly starting from the flat configuration and 
measuring xl a plateau at the value of M(/3) is 
reached. To be sure that the initial topology is 
not changed by the heating procedure, each con- 
figuration of the sample is checked during heating 
by performing a few cooling steps to detect its 
topological charge on a copy of it. If a change 
occurred the configuration is discarded [p^ . 

Once M and Z are known, x can be extracted 
by a subtraction procedure following eq. (|^). 
Clearly Z and M depend on the operator used 
for Ql- a standard definition for is |l9j 

-1 ^'^ 

QUx)^^^ e^^p^Tr (n^^Hp^) . (8) 



Here e^j/po- is the standard Levi-Civita tensor for 
positive directions while for negative ones the re- 
lation ifj^^pa- = -(■-fj.^pa holds. is the pla- 
quette in the fi — v plane at point x. With this 
definition Z 2± 0.18 and the mixing M is most 
of the measured signal in the scaling region [Tl] ], 
so that a large error derives from the subtraction 
procedure performed to extract x- 

In order to obtain a more precise determina- 
tion of X, a new definition of can be used for 
which M is strongly reduced and Z is closer to 
its continuum value Z — 1. 

In this work we use the sequence of improved 
operators Q^*'' defined in ref. the starting 

operator of the sequence is the standard defini- 
tion of eq. (H) and a smearing procedure is per- 
formed at each step of the sequence. We stop 
at the second smearing, where a considerable im- 
provement is already achieved. For a full account 
of our results see ref. . 

3. RESULTS 

We measured the topological susceptibility at 
zero temperature on a symmetric lattice 16'* 
and at finite temperature on a lattice 32"^ x 
8. The simulations were performed on a APE 
QUADRICS machine. 

In figure 1 we plot the value of (x)^^/*-* at zero 

(i) 

temperature versus /3 for the three definitions Xl 
{i = 0, 1, 2) of the lattice topological susceptibil- 
ity. There is good scaling: {xY"^^^^ is independ- 
ent of f3, as it should. There is also an excellent 
agreement between the three determinations. To 
fix the value of (x)*-^/*^ in physical units we use 
for a{f3) the two loop formula and for A^ the de- 
termination of ref. HI] Al = 4.56(11) MeV. The 
horizontal line is the linear fit to the 2-smear data. 
It gives (x)^^^'*^ = 175(5) MeV and is consistent, 
within errors, with that of ref. JTH-p^. The error 
includes the uncertainty in A^. 

In figure 2 the topological susceptibility x = 
(Xl' - M(*))/(Z(*)2^4) ^Yie transition point is 
shown for the 1 and 2-smeared operator, x drops 
by one order of magnitude from the confined to 
the deconfined phase. The results obtained with 
the two operators are compatible as they should. 
The data for the 0-smeared operator have very 
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Figure 1. x at T = 0. The straight hne is the 
result of the hnear fit of the 2-smeared data. The 
improvement from Q^^^ to Q^^^ is clearly visible. 
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Figure 2. x/A|^ versus T /Tc across the deconfin- 
ing phase transition. The horizontal band is the 
determination at T = of figure 1. 



large error bars and are not shown in the figure. 
The data have been plotted versus T/Tc where 
Tc is the deconfining temperature coresponding 
to Pc{Nr = 8) = 6.0609(9) @]. To determine 
T/Tc we only need the ratio a{(3c)/a{(3) and for 
that the two-loop expression is certainly a good 
approximation within the small interval of /3 used, 
where A^, can be considered as a constant. 

The solid line in figure 5 corresponds to the 
value of X zero-temperature and is consistent 
with the data below Tc indicating that x is prac- 
tically T-independent in the confined phase. 

We thank Graham Boyd and Enrico Meggi- 
olaro for useful discussions. A.D.G. acknow- 
ledges an interesting discussion with Heinrich 
Leutwyler. 
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